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ABSTRACT 

We used near-infrared 2MASS data to construct visual extinction maps of a 
sample of Southern Bok globules utilizing the NICE method. We derived radial 
extinction profiles of dense cores identified in the globules and analyzed their 
stability against gravitational collapse with isothermal Bonnor-Ebert spheres. 
The frequency distribution of the stability parameter {^max) of these cores shows 
that a large number of them are located in stable states, followed by an abrupt 
decrease of cores in unstable states. This decrease is steeper for globules with 
associated IRAS point sources than for starless globules. Moreover, globules in 
stable states have a Bonnor-Ebert temperature of T = 15 ± 6 K, while the group 
of critical plus unstable globules has a different temperature of T = 10 ± 3 K. 
Distances were estimated to all the globules studied in this work and the spectral 
class of the IRAS sources was calculated. No variations were found in the stability 
parameters of the cores and the spectral class of their associated IRAS sources. 
On the basis of ^^CO J = 1-0 molecular line observations, we identified and 
modeled a blue-assymetric line profile toward a globule of the sample, obtaining 
an upper limit infall speed of 0.25 km s~^. 

Subject headings: dust, extinction — infrared: ISM — ISM: globules — stars: 
formation 
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Introduction 



Bok &: Reillyl (119471 ) were the first to call attention to small dark regions, dense and 
rounded, which called globules, and suggested the possibility that the star formation could 
occur within them. These objects were called later Bok globules and in general they are 
isolated, cold (T ~ 10 K), dense (n ~ 10^~^ c ni~^), and seen as pote ntial sites for the 
formation of low-mass stars (< 10 M©). Following iBeichman et al.l (119861 ). the globules that 
do not contain protostars or young stellar objects (YSOs) associated, are called starless cores. 
When they show s ome evidence of st ar formation, such as Class or I protostars, are called 
protostellar-cores (I Andre et al.lll993l ). 



As the majority of the stars form in molecular cloud complexes and in clusters embedded 
in giant molecular clouds and not in isolation, it is very important to know in detail the 
physical processes that give rise to the formation of a single star in an isolated environment. 
For this purpos e, Bok globules const itute an ideal laboratory to study the isolated low-mass 
star formation (IBourke et al.lll995al . hereafter BHR). 



Visual extinction was used by lClemens fc Barvainid (119881 ) to study 248 Bok globules in 
the Northern Hemisphere, and by BHR to study 169 globules in the Southern Hemisphere. 



Sub-m illimeter and millimeter wavelength studies of starless cores (jWard-Thompson et al. 



19941 ). indicate that the ID profiles (mean of azimuthal distributions) are flat for radii ~ 
3000-7000 AU, and steeper for larger radii. These profi les d e viate from the distribution 
n(r) oc of the singular isothermal sphere, proposed by IShul (119771 ) as the initial state of 
isolated dense cores before gravitational collapse. Similar profiles are obtained by mapping 
the dust distribution in a den se core or Bok globule. T hese configurations are consistent 
with "Bonnor-Ebert spheres" (jEbertlll955l : lBonnorlll956l ). which are non-singular solutions 
to the equations of hydrostatic equilibrium confined by an external pressure. In the ID 
isothermal approximation, there is a family of Bonnor-Ebert solutions parametrized by the 
central density Uc, where the density profile is characterized by two regimes: a slow decrease 
in the density for small radii, and a faster decrease, of power-law type (oc r"^), for larger 
radii. 



Detailed modelling of pr ofiles in the infrared (jKandori et al.ll2005l : iTeixeira et al.ll2005l ) 



and continuum emission (e.g.. lWard-Thompson et al.lll999l ) confirm that Bonnor-Ebert spheres 



are a good approach to the internal structure of starless cores. However, it is not currently 
possibly to distinguish between a static Bonnor-Ebert sphere and a collapsing one for profiles 
based in dust absorption or emission, since the structur e of the profi le, for a given ric, does 
not change significantly until the last stages of collapse (lMyersll2005l ). 



Using near-infrared observations, lAlves et al.l (120011 ) showed that a Bonnor-Ebert sphere 
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with stability parameter ^rnax = 6.9 adjust pe r fectly the observed column density profile of 
the starless globule Barnard 68. iHarvev et al.l (l200l[) ra odeled the profile of the protostellar 
globule B335 with ^^ax = 12.5, and lHarvey et al.l (I20031) estimated ^mn-r = 25 for the starless 
core L694-2. The Coalsa c k glob ule 2 was modeled by iRacca et al.l (120021 ) that derived ^rnax 
7.3, while iLada et al.l (120041) o btained ^max = 5.8 using a more sensitive telescope. In 



Kainulainen et al 



7 _ ^ ^ * < (IJyii <3J - - 

recent works. iHuard et al.l (120061 ) found ^max = 35.8 for the protostellar core L1014, and 



(l2007h obtained ^„ 



23 for the protostellar globule DCld303.8-14.2 



and ^rnax ^ 8 for the starless globule Thumbprint Nebula. 



The largest sample of Bok globules studied in the near-infrared was done by lKandori et al 



( 120051 ) where 14 globules were selected. They found that more than half of the starless glob- 
ules (7 out of 11) a re located close to the critical state, with ^^nax = 6.5 ± 2. This led 
Kandori et al. (j2005 ) to suggest that a Bonnor-Ebert sphere in the critical state character- 
izes the structure of typical starless cores. The remaining starless globules and those that 
present evidence of star formation show clearly unstable states, with C,max > 10. 

In this work, we present a study of 21 Southern Bok globules, where 11 are starless and 
10 have associated IRA^fl point sources. In §2]we present the source list and the observations. 
Using 2MASS data, we constructed visual extinction maps in ^ which enabled us to detect 
dense cores embedded in the globules. In §1] we estimated new distances to the globules, 
and by means of the Bonnor-Ebert model, in §5] we studied their internal structure and 
stability against gravitational collapse. In ^ we analized the IRAS point sources associated 
to the globules and in §7] we found evidence of infall motions in a globule using millimetric 
observations of the ^^CO and C^^O molecules. Our main conclusions are summarized in ^ 



2. Observations 
2.1. Sample List 

In this work, we searched for a list of Bok globules satisfying the following conditions: 

1. that they were completely isolated, i.e., neither associated to bright nebula nor to a 
molecular complex. 



2. that they had a distance determination. 



^Infrared Astronomical Satellite (jNeugebauer et al.lll984l ) 
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Using these criteria, we got 21 Southern Bok globules, selected among 169 globules of the 
BHR sample, where 11 do not present any evidence of star formation, and 10 have associated 
IRAS point sources. Following previous definitions (see in the first group there are 
"starless globules" , and in the second there are "IRAS globules" . 

Table [1] presents the sel ected objects. Column 1 and 2 show the names of the globules 
as they appear in BHR and [Hartley et al.l (119861 ). respectively; columns 3 and 4 give their 
equatorial coordinates; column 5 shows the optical sizes in arcmins, and column 6 gives the 
distances determined by BHR. 



2.2. Near-Infrared Data 



Th e near-infrared data w ere obtained from the Two-Micron All Sky Survey (2MASS) 
catalog ( Skrutskie et al.|[2006l ). The 2MASS Point Source Catalog (PSC) contains 471 million 



stars, covering 99.998 % of the sky in the J(1.24 yum), H(1.66 /im) and Ks(2.16 /xmj^ bands 
observed with two telescopes of 1.3 meters diameter, located in Mount Hopkins (Arizona, 
USA) and Cerro Tololo (Chile). The 2MASS PSC is more than 99% complete for J < 15.8, 
H < 15.1 and K < 14.3. 

For each of the globules in Table [H we retrieved JHK magnitudes and their correspond- 
ing uncertainties in the 2MASS PSC in regions of 15' x 15' centered on the globules. Only 
those stars whose photometric uncertainties were < 0.1 mag in all three bands were ex- 
tracted, ensuring a signal-to-noise ratio S/N > 10. The search for stars in the 2MASS PSC 
was made through the web interface GatoiE. 



2.3. Molecular Line Data 

We searched for molecular lines observed toward these globules having line profiles data 
available. We found only ^^CO and C^^O (J = 1-0) line observati ons toward the globule s 



BHR 138 and BHR 149 made in 1992 October and 1993 May fsee IVilas-Boas et al.lll994[ ). 
These observations were made with the 15 m Swedish-ESO Submillimeter Telescope (SEST) 
at La Silla, Chile. The receiver front end was based on a Schottky diode waveguide mixer 
followed by an intermediate-frequency amplifier. The SSB system temperature was typically 



^From hereafter, filter "Kg" will be called simply "K". 



3 



http : //irsa. ipac . caltech.edu/applications/Gator/ 
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Table 1. Sample List 



Name DC No. R.A. Dec. Size Distance 

(J2000) (J2000) (a' x 6') (pc) 



Starless Globules 



016 


255.4-3.9 


08 05 26 


-39 08 54 


6 


X 


2 


300 


044 


269.5+4.0 


09 26 19 


-45 11 00 


2 


X 


2 


300 


053 


274.2-0.4 


09 28 47 


-51 36 42 


8 


X 


4 


500 


059 


291.1-1.7 


11 07 07 


-62 05 48 


5 


X 


3 


250 


074 


300.0-3.7 


12 22 09 


-66 27 06 


3 


X 


3 


175 


075 


300.2-3.5 


12 24 13 


-66 10 42 


6 


X 


3 


175 


111 


327.2-1-1.8 


15 42 20 


-52 49 06 


3 


X 


2 


250 


113 


331.0-0.7 


16 12 43 


-52 15 36 


8 


X 


4 


200 


133 


340.5+0.5 


16 46 45 


-44 30 48 


10 


X 


4 


700 


144 


346.4+7.9 


16 37 28 


-35 13 54 


7 


X 


1 


170 


145 


347.5-8.0 


17 48 01 


-43 43 12 


8 


X 


5 


450 


IRAS Globules 


034 


267.2-7.2 


08 26 34 


-50 39 54 


8 


X 


2 


400 


058 


289.3-2.8 


10 49 00 


-62 23 06 


5 


X 


3 


250 


117 


334.6+4.6 


16 06 18 


-45 55 18 


7 


X 


2 


250 


121 


337.1-4.9 


16 58 42 


-50 35 48 


2 


X 


2 


300 


126 


338.6+9.5 


16 04 29 


-39 37 48 


6 


X 


5 


170 


138 


345.0-3.5 


17 19 36 


-43 27 06 


10 


X 


6 


400 


139 


345.2-3.6 


17 20 45 


-43 20 30 


4 


X 


2 


400 


140 


345.4-4.0 


17 22 55 


-43 22 36 


6 


X 


3 


400 


148 


349.0+3.0 


17 04 26 


-36 18 48 


2.5 


X 


2.5 


200 


149 


349.2+3.1 


17 04 27 


-36 08 24 


3 


X 


1 


200 



Note. — Units of right ascension are hours, minutes, and seconds, 
and units of declination are degrees, arcminutes, and arcseconds. 
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390 K. 

The back end used an acousto-optical spectrometer with resolution 43 kHz (0.11 km 
s~^) and total bandwidth of 100 MHz. The spectra were taken by using overlap frequency 
switching mode (with 7 MHz frequency shift) and integrating, on average, for periods of 2 
minutes. The observations were chopped against a cold load to obtain the correction for 
atmospheric attenuation. The halfpower beam width was 48" and the beam efficiency 0.9. 
The rms pointing accuracy was better than 10". 



3. Visual Extinction Maps 



We used the NICE method of iLada et al.l (jl994j ) to construct visual extinction maps 
for the globules. This technique uses a control field region free of reddening to refer the 
extinctions measured in the cloud region, and assumes that the stellar population is the 
same in the control field and the cloud region. If we know the intrinsic color of a star, then 
the color excess is: 



E{H -K) = {H -K) 



observed 



-(H-K) 



intrinsic- 



Intrinsic {H — K) colors for main-sequence and gian t stars are in a s r nall interval, 0.0 to 0. 3 
for stars with spectral types between AO and M5-7 (iKoornneeil 1 1 9831 : iBessell fc Brettlll988l ). 
Therefore, we can approximate the intrinsic color of each on-cloud star by the average color 
of the control field, and the color excess is: 



E{H - K) = {H - K) - {H - K) 



control ' 



(2) 



Using equation [21 we calculate the color excess for each star in the globule region. 
Then, we can express this color excess in terms o f visual extinction with the knowledge of 
the reddening law in the 2MASS system. Folowing lNielbock fc Chinil (120051 ). this relation is: 



Av = l9Ax E{H - K), (3) 

whose associated uncertainty is: 



aAy = 19.4 X J erf + al^t.^i, (4) 
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where (Tj is the uncertainty in the observed [H—K) color and (Tcontroi is the standard deviation 
of the colors in the control field. Nearly 90% of the regions have (Jcontroi ~ 0.1 mag, which 
sets the minimum mensurable value for the extinction in ~ 2 mag. Only two regions in our 
sample have (Jcontroi ~ 0.2 mag, giving a minimum of ~ 4 mag of visual extinction. 

In this way, we obtained the distribution of extinction in the fields containing the Bok 
globules. However, the spatial distribution of stars is non-uniform. In order to produce a uni- 
formly sampled extinction map, we smoothed the data by spatially convolving the extinction 
measurements (equation [3]) with a Gaussian filter (kernel) with a given sm oothing paramete r 



(resolution), and finally sampled the map at the Nyquist frequency (e.g.. lLada et al.lll999l ). 
The form of the Gaussian kernel is given by: 



K{a,ai,6,6.i) = —expi-^] , (5) 



where 



r^^ = (5_5,)^ + (a-a,)W(5) 



(6) 



and h is the map resolution. Hence, the visual extinction Ay at each point {a, 6) in the map 
is calculated using: 



Av{a, 6) 



(7) 



where n is the total number of stars observed in the cloud region and Ay(aj, 5i) is given by 
eq uation [3l and the associa ted uncertainty for the extinction at the point (a, 5) in the map 
is dlombardi fc AlveslboOlh : 



^Er=i -^^(Q'^ Si) X a\^iai, 5i) 



where (T^^(aj, 6i) is given by equation HI 



3.1. Construction of the Extinction Maps 



To test the 2MASS data using th e NICE met h od, w e applied it to Barnard 68 (B68), a 
well known and studied Bok globule. lAlves et al.l (l200ll ) observed B68 with the ESO-NTT 
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telescope in the JHK bands to produce an extinction map of this globule with a resolution 
of 10". 

The choice of the smoothing parameter h is a compromise between a final map with 
high signal-to-noise ratio and low spatial resolution, or a noisy map with high resolution 
( iLombardi fc Alvesll200ll ). To construct the extinction map of B68 (and all the globules in 
this work) we choose h = 20". Nonetheless, other values were tested: with h = 10" the 
map is very noisy, and with h = 30" the map loses details. Therefore, a choice oi h = 
20" provides the best spatial resolution for a signal-to-noise ratio between 8 and 12 in the 
central region of B68. Figure [1] shows the visual extinction m ap for B68 constru cted using 
data from the 2MASS catalog. Comparing with the map of lAlves et al.l (120011 ) . we see a 
good correspon dence betwee n them, taking into account that the 2MASS used a 1.3 meters 
telescope, while lAlves et al.l used a 3.5 meters telescope and longer integration times. The 
sensitivity difference between these telescopes is certainly responsible for the different values 
of extinction found in the center o f B68. From the 2MA SS data, we obtain a maximum 
extinction of 24 magnitudes, while lAlves et al.l (120011 ) and iHotzel et al.l (l2002bl ). using the 
same data set, derived 30 magnitudes. This result shows that the extinctions derived toward 
the center of the globules in our sample, as well as the density contrast, have to be seen as 
lower limits. In despite of getting a lower extinction peak, the extinction distribution and 
the peak position are coincident. 

Figure [2] shows the visual extinction maps for all Bok globules in our sample. They 
were constructed with the NICE method as described in ^ Each map has ~ 15' x 15', 
centered on the position listed in Table [H In the right side of each map, a bar indicates the 
extinction in magnitudes. White crosses superimposed denote the IRAS sources associated 
to the globules (see 



3.2. Comparison with the Star Counts Method 

In this item we compare star counts and the NICE methods in order to verify which 
one is more apropriate to study extinction toward the Bok globules of our sample, derived 
from 2MASS data. 

The 2MASS maps constructed in this work were performed only with those stars whose 
color excess E{H — K) > 0, thus removing almost all foreground stars that can be present 
in each globule region. This allowed us to minimize the effects th at foreground stars ca n 



produce in the construction of the extinction maps (see Figure 6 of ICambresy et al.ll2002l ) 



In order to test if the 2MASS data together with the NICE method can produce reliable 
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dust extinction maps of Bok globules, we also constructed maps using the traditional star 
counts method. The number of stars per unit area brigher than rrix in the control field, Nofj, 
is: 



\og{Noff) = a + bmx. (9) 

On the cloud region, the slope of equation [H] remains the same, but the number of stars per 
unit area, iVo„, decreases as extinction increases. The extinction is then: 

A. = hos(^), (10) 



b "\ N, 



wher e b is the slope of the cumulative luminosity function (see, e.g., iGregorio Hetem et al. 



19881 ). We converted the Ax extinctio n (with A = J,H,K) to o ptical Ay extinction using the 



reddening law in the 2MASS system (INielbock &: Chinill2005l ). 



We constructed extinction maps for B68, our test globule, and for BHR 059, one of the 
smallest globules in our sample, using the star counts technique with the condition that, at 
least, one star must be present in each counting box. We obtained three Nyquist sampled 
maps for each globule, that is, one map for counts in J, H and K bands, respectively. Table 
|2] shows the parameters used in equation [10] to build the maps. 

The extinction maps for B68 and BHR 059 are shown in Figure [3l We can see that 
B68 is not detected in the map based on K band star counts. This is probably because the 
distribution of the stars in the K band is almost uniform in the most opaque region of B68. 
The H and J band maps show B68 as a region of extent ~ 3-4 pixels without any details on 

Table 2. Parameters of the Star Counts Maps 



Barnard 68 BHR 059 



Wavelenght Bin Nofj b Bin Nojf b 



J 


80" 


21 


0.36 ± 0.03 


70" 


16 


0.35 ± 0.02 


H 


50" 


8 


0.35 ± 0.03 


70" 


15 


0.34 ± 0.01 


K 


50" 


7 


0.37 ± 0.02 


70" 


13 


0.32 ± 0.01 
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the distribution of dust extinction in the core. It is important to note that a non uniforme 
distribution of field stars causes unreahstic peaks of high visual extinction, which are not 
seen in the DSS red plates. The same situation occurs for BHR 059. 



We estimated the erros in the derived visual extinctions (jDickmanlll978l ) and obtained 
SAy ~ 4, 8 and 12 magnitudes for J, H and K band maps, respectively, in regions where 
Ay = (where Non = ^o//) and 6 Ay ~ 10, 20 and 30 magnitudes for J, H and K band 
maps, respectively, in the most opaque regions (where Non = 1 in a single counting box). 
On the other hand, the signal-to-noise ratio for the maps constructed through the NICE 
method is always greater than 3 in the most extincted regions where we constructed the 
radial extinction profiles. Even though the signal-to-noise ratio is not as good for the more 
diffuse, low extinction regions, it is not crucial for our purposes, because we are interested 
in the profile of the core embedded in the globule to finally model it with a Bonnor-Ebert 
sphere. This discussion demonstrate that, to our purposes, the NICE method applied to the 
2MASS data is the best way to construct visual extinction maps and trace Bok globules to 
finally apply the stability analysis of the radial profile. 



4. Distance Determination 



The knowledge of the distances to the globules is extremely important, as they are 
necessary to determine several fundamental parameters, such as mass, size, density and 
also to calculate the luminosities of the young objects embedded in the globules. Several 
methods have been used to determine distances , being th e more typical the star count 
me thod (IBok &: Bokll941l) and the Wolf diagram (jWolilll923l ). A newer technique, proposed 
bv lMaheswar et all J20o'l . uses optical and near-infrared broadband photometry of stars in 
the field of the cloud. 



In this paper, we employ the technique of iDickman fc Clemens! (Il983l ). which uses the 
plot of color excess E{B — V) versus the distance to stars in the vicinity of the cloud. In 
this graph, a jump in the reddening indicates the presence of a cloud or a group of clouds. 
The same method was used by BHR to estimate the distance to some globules. However, 
for the reasons described in §4.11 we decided to calculate ourselves these distances again. 



4.1. Application of the Method 



To estimate the distances to the Bok globules , we need to know the redde ning E{B — V) 
of the stars in a given area around each globule. iDickman fc ClemensI ( 119831 ) used a region 
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of 11° X 16° in their work. BHR used circular regions of 5° radius. If the region had not 
a sufficient amount of stars, they increa s ed th e radius to 7.5° or 10°. In this work, we 
used the SKY2000 catalog of iMyers et al.l (|2002| ). which has information for ~ 300,000 stars 
brighter than 8 magnitudes, and searched for stars in circles of 3° radius centered on each 
Bok globule. Figure H] shows the galactic distribution of the globules, and the circles are 
their neighborhood used to determine the distance. As we see, some regions coincide, given 
the proximity between the globules. In these cases, both are assumed to be at the same 
distance. 

To calculate the color excess E{B — V) for each star, we searched the SKY2000 catalog 
for all those stars with known spectral types, apparent visual magnitude V, and observed 
color index {B — V). Having the spectral type of a star, we know t he absolute ma g nitud e 
My and the intrinsic color index [B — V)q using the calibration of ISchmidt-Kalerl ( 1982 ). 



Assuming a value of 3.1 for the total-to-selective extinction ratio (IRieke fc Lebofsky 



19851) 



we calculate the extinction along the line of sight of a star in the neighborhood of a globule 
as: 



Ay = 3.1 X E{B -V) = 3.1 X [{B - V) - {B - V)o] , 



and finally the distance in parsecs, corrected for extinction: 



log(r) 



-{V - My + 5- A 



V 



(12) 



Figure \5\ shows the graphics of E{B — V) versus distance for the globules, where the 
dashed lines indicate the distances adopted here. In Table [3l column 1 indicates the BHR 
globule name, column 2 the distance obtained in this work, and column 3 the distance 
obtained by BHR. 



5. Internal Structure: Bonnor-Ebert Models 



To study the internal structure of the Bok globules, they were fitted with Bonnor-Ebert 
isothermal spheres. From the extinction maps generated in ^ radial extinction profiles were 
built for all the Bok globules considered here. 

A sphere of self-graviting gas in hydrostatic equilibrium, where the pressure and density 
at each point are related through t he isotherrnal equ ation of state, is called a Bonnor-Ebert 
sphere ( lEbertl Il955l : iBonnorl Il956l ). iBonnorl (119561 ) studied the gravitational stability of 



- 12 - 



Table 3. Distances to the Bok Globules 



Name 


Distance^ 


Distance'^ 




(pc) 


(pc) 


016 


250 


300 


no A 


zOO 


A AA 

400 


U44 


200 


OA A 

300 




200 


C AA 

500 


058 


200 


250 


059 


O AA 

200 


OCA 

250 


074 


175 


1 TC 

175 


075 


175 


175 


T T 1 

111 


OCA 

250 


o cr A 

250 


113 


200 


200 


117 


175 


250 


121 


125 


300 


126 


250 


170 


133 


225 


700 


138 


225 


400 


139 


225 


400 


140 


225 


400 


144 


225 


170 


145 


150 


450 


148 


175 


200 


149 


175 


200 


'^Distance obtained in 


this work. 


''Distance 


obtained 


bv Bourke et al. 



(11995ar i. 
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isothermal spheres and deducted a critical parameter that allows classify them as stable or 
unstable against gravitational collapse. 

The shape of a Bonnor-Ebert density profile consists of a fiat central region and a steeper 
external region, of the form p(r) oc r^^. For a Bonnor-Ebert sphere with radius R and central 
density pc, its normalized profile, i.e., p{r) / pc vs. r/i?, is characterized by a dimensionless 
radius ^max, whose critical value is 6.5. This value of the critical parameter corresponds 
to a contrast between the central density and the density at the edge of the sphere of 14. 
Increasing values of this parameter denote spheres more centrally condensed. 



5.1. Radial Extinction Profiles, Masses and Densities 

In order to construct radial profiles, we define the position of the center of the globule 
as the position of the extinction peak, which does not necessarily coincide with the position 
listed in Table [H As some extinction maps have more than one core for the same globule, 
we will assign a letter to the new cores identified within a single globule. Then, we set 
concentric annuli at each core center, 20" wide, and averaged the extinction at the Nyquist 
frequency, obtaining a mean extinction value every 10" (impact parameters). The resulting 
uncertainty for an impact parameter corresponds to the error propagation of each pixel 
within the annulus, calculated with equation [HI The radial profile for B68 is shown in Figure 
[6l and the resulting radial extinction profiles for the globules of our sample are shown in 
Figure [71 The solid curve in each plot corresponds to the Bonnor-Ebert fitting (see §5.2p . 

As shown in §5.21 and §5.31 below, to fit these profiles we modeled a radial density profile 
for each impact parameter, and the radiu s R of the core i s the impact parameter for which 



we obtained the best fit for that profile. iTeixeira et al.l (120051 ) define the radius of a core 



as the distance for which the profile reaches a constant lev el of extinction or plateau. The 



radius obtained in this article coincide with the definition of ITeixeira et al.l because the cores 
are embedded in a more diffuse region in the globules. Table [H lists in column 1 the BHR 
globules, where a letter attached represents a substructure identified in the globule, columns 
2 and 3 give the coordinates of the cores as defined above, and column 4 lists their sizes. 

Once constructed the radial profiles of the cores, we calculate their masses integrating 
the profile until the radius of the core: 



M = pniH / Nh dn, (13) 
Jn 

where /i (= 1.36) is the mean atomic weight of the gas, rriH is the mass of the H atom, is 
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Table 4. Coordinates and Sizes of Dense Cores in Bok Globules 



Name 


R.A. 


Dec. 




R 




(J2000) 


(J2000) 


(pc) 






Slaiicss 


Globule 








016 A 


08 05 18.8 


-39 


08 


53 


0.07 


016 B 


08 05 27.2 


-39 


08 


20 


0.06 


044 A 


09 26 09.3 


-45 


11 


10 


0.07 


044 B 


09 26 20.4 


-45 


10 


56 


0.06 


053 


09 28 46.5 


-51 


36 


20 


0.07 


059 


11 07 10.1 


-62 


05 


36 


0.06 


074 


12 22 10.2 


-66 


27 


39 


0.07 


075 


12 24 14.2 


-66 


10 


58 


0.09 


111 


15 42 19.8 


-52 


48 


26 


0.11 


113 


16 12 51.6 


-52 


16 


23 


0.07 


133 


16 46 42.6 


-44 


31 


10 


0.13 


144 A 


16 37 29.3 


-35 


13 


43 


0.07 


144 B 


16 37 35.5 


-35 


14 


43 


0.07 


145 


17 47 51.9 


-43 


42 


15 


0.07 


IRAS Globules 


034 


08 26 27.7 


-50 


39 


30 


0.06 


058 


10 49 02.4 


-62 


22 


18 


0.05 


117 A 


16 06 25.2 


-45 


54 


17 


0.06 


117 B 


16 06 11.8 


-45 


56 


29 


0.06 


121 


16 58 47.0 


-50 


36 


34 


0.06 


126 


16 04 29.2 


-39 


37 


47 


0.11 


138 


17 19 32.9 


-43 


26 


55 


0.10 


139 


17 20 51.7 


-43 


19 


44 


0.07 


140 A 


17 22 56.2 


-43 


22 


26 


0.07 


140 B 


17 22 52.9 


-43 


21 


40 


0.06 


140 C 


17 23 29.6 


-43 


25 


10 


0.06 


148 


17 04 26.4 


-36 


18 


34 


0.07 


149 A 


17 04 31.2 


-36 


07 


52 


0.06 


1 i<) I! 


17 01 5:-!.! 


-:5() 


0:-! 


15 


0.0() 



Note. — Units of right ascension arc 
hours, minutes, and seconds, and units of 
declination are degrees, arcminutes, and 
arcseconds. 
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the column density and is the cloud area projected in the plane of the sky. T o express the 



column density in terms of visual extinction, we used the gas-to-dust ratio of iBohlin et al. 



(119781): 



Nh = 2x 10^^ Av cm-2 mag-\ (14) 

and substracting the contribution of the diffuse structure or plateau extinction, we have that 
the mass of an embedded dense core is: 

M = 2n (^^^ fimnD^ £ (^Ay - A^"*^"") rdr, (15) 

where D is the distance to the globule and R is the core radius. Finally, to calculate the 
mean volumetric density of each core, we used: 

3M , , 

The above results are listed in Table [51 Column 1 indicates the core name, column 
2 the plateau extinction, column 3 and 4 are the masses and mean volumetric densities of 
the cores without plateau substraction, and in columns 5 and 6 are given the masses and 
densities substacting the plateau contribution. 



5.2. Bonnor-Ebert Model 

A Bonnor-Ebert sphere is a pressure-truncated, self-graviting isothermal sphere of gas in 
hydrostatic equilibrium. The equation that describes such a gaseous sphere is the modified 
Lane-Emden equation: 

where p = pce~'^^^\ Here, ^ = {r / a)^ AixGpc is a dimensionless radial parameter and a 
is the isothermal sound speed. Introducing standard boundary conditions (0(0) = and 
d0(O)/d^ = 0), we solved equation 1T71 numerically using a fourth order Runge-Kutta method. 



Assum ing an external radius R, solutions of equation 1T71 can be parametrized by pc- iBonnor 



( Il956l ) demonstrated that the value of ^ evaluated at the external radius R: 
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Table 5. Masses and Densities of Dense Cores in Bok Globules 



Name 














(mag) 


{Mq) 


(10* cm-3) 




(10* cm-3) 




Starless Globules 


016 A 


8.7 


9.5 


19.1 


3.6 


7.2 


016 B 


6.2 


4.7 


15.1 


1.2 


3.8 


044 A 


5.1 


5.4 


10.9 


1.4 


2.7 


044 B 


5.8 


4.0 


12.9 


1.1 


3.4 


053 


7.1 


8.0 


16.1 


2.4 


4.9 


059 


8.1 


6.8 


21.7 


2.6 


8.4 


074 


2.8 


3.6 


7.4 


1.5 


3.0 


075 


3.5 


6.3 


5.9 


1.8 


1.7 


111 


9.8 


22.7 


11.8 


5.5 


2.9 


113 


7.9 


7.7 


15.5 


1.4 


2.9 


133 


7.5 


20.4 


6.4 


3.2 


1.0 


144 A 


11.7 


11.5 


23.2 


2.6 


5.3 


144 B 


11.1 


9.4 


18.9 


2.2 


4.4 


145 


3.6 


3.9 


7.8 


1.4 


2.8 


IRAS Globules 


034 


2.5 


2.2 


7.0 


0.9 


2.8 


058 


5.4 


4.4 


24.3 


2.1 


11.4 


117 A 


8.9 


4.9 


15.8 


0.9 


2.8 


117 B 


8.3 


5.9 


18.8 


0.8 


2.7 


121 


3.1 


2.4 


7.8 


0.7 


2.4 


126 


6.0 


13.0 


6.8 


2.6 


1.3 


138 


1.9 


5.0 


3.5 


2.3 


1.6 


139 


4.8 


6.3 


12.8 


2.7 


5.4 


140 A 


7.9 


6.8 


13.7 


1.6 


3.3 


140 B 


7.7 


5.2 


16.6 


0.9 


3.0 


140 C 


4.5 


3.8 


12.3 


1.4 


4.4 


148 


4.5 


4.8 


9.8 


1.4 


2.8 


149 A 


4.6 


4.0 


12.8 


1.5 


4.9 


149 B 


3.6 


2.8 


8.9 


0.9 


2.8 



''Value without plateau substraction. 
''Value with plateau substraction. 



provides a stability measure. Systems with ^max > 6.5 are unstable to gravitational collapse. 
Equivalently, the density contrast between the center and the border of the cloud is a function 

of ^max- 




(19) 



PR 

If pr = p{R), and in the critical state C,max = 6.5, then pd Pr = 14. Equation [T^ shows that, 
for stable states where ^max < 6.5, the center-to-edge density contrast is lesser than 14. 

Finally, knowing the core radius R and the stability parameter ^max-, we can calculate 
the physical parameters of the sphere, as the central density: 



or the volumetric central density: 



where p (= 2.33) is the molecular weight of the H2 molecule. The mass of the Bonnor-Ebert 
sphere is given by: 




(20) 




(22) 



and the external pressure at the edge: 



ext 




(23) 



5.3. Bonnor-Ebert Fitting 



In order to compare theoretical profiles with observed extinction profiles for each core, 
we constructed a series of theoretical Bonnor-Ebert profiles in the following way: for various 
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values of the stability parameter ^max and temperature T, we obtained the volumetric density 
n{r) solving equation [T7| as explained above. Then, we integrated the volumetric density 
profile along the line-of-sight to obtain the hydrogen column density profile: 



NsEir) = 2 X 



^ , r'dr' 
n[r j 

V — ; 



(24) 



where r is the projected distance from the center of the core or impact parameter. The 
stability parameter was varied from 3 to 15 in steps of A,^ = 0.1, and the temperature was 
varied from 5 K to 30 K in steps of AT = 0.1 K. To expres s the column density in terms 
of visual extinction, we used again the gas-to-dust ratio of iBohlin et al.l (119781 ) given by 
equation [HI The fitting was evaluated using a reduced with 2 degrees of freedom: 



n 



E 

i=l 



A 



BE/ 



H 2 



0",; 



(25) 



where n is the number of points considered in the profile and is the uncertainty in the 
observed Ay. To estimate the uncertainty of the best-fit parameters, we took the second 
derivative of Xr with respect to each parameter ^max and T in the region where Xr is minimum: 



^max \ pi Cl 



(26a) 



2 ^ {d'x^ 
arp = z 



according to iBevington fc RobinsonI (119921 ) 



i2,,2\ -1 



9T2 



(26b) 



The Bonnor-Ebert model applied to B68 using the 2MASS data gave the following 
results: 



'imax — 6.9 ih 0.3, 



T = 16.8 ± 0.9 K. 



The theoretical profile for B68 is shown in Figure as a continuous line superimposed on 
the observed profi le. The stability parameter is in good concordance with that obtainded by 
Alves et al.l (I20011). w ho calculated ^max = 6.9 ±0.2 using a temperature of 16 K measured by 
Bourke et al.l (Il995bl ). which coincides with the value determined in our fit. However, more 
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recent observations indicate a lower value for the temperature of B68. iHotzel et a 



observ ed B68 in NH3 and obtained a gas kinetic temp erature of 10 ± 1-2 K , while 



(2002a) 



Lai et al. 



( 120031 ) found a value of 11 K using the same molecule. iHotzel et al.l (l2002bl ). using the same 
data, modeled B68 with a Bonnor-Ebert sphere and found central extinctions of Ay = 30.3 
mag, which is different from the value derived in this work. Ay = 19.5 mag. This difference 
arises probabely because they observed B68 with a greater resolution and sensitivity than 
the 2MASS, making it possible to detect a larger number of stars toward B68 than 2MASS. 

Table [6] summarizes the results of the Bonnor-Ebert modelling for the Bok globules of 
our sample. Column 1 gives the names of the cores, column 2 the stability parameters, 
column 3 the Bonnor-Ebert temperatures, column 4 and 5 the central volumetric densities 
and masses, column 6 the external pressures in units of the Boltzmann constant, column 7 
the central extinctions and column 8 the minimum value of Xr- The resulting theoretical 
profiles for the globules are shown in Figure [71 

For the physical parameters obtained in the Bonnor-Ebert fitting, we calculated the 
mean, standard deviation (SD), median, lower quartile (LQ), and upper quartile (UQ) for 
the whole sample, for the starless globules and for the IRAS globules, respectively. These 
quantities are shown in Table [3 Comparing starless with IRAS globules, we see from this 
table that the mean value of the physical parameters are nearly identical for both kind of 
globules. However, for the temperature and external pressure, starless globules have larger 
values than IRAS globules. In order to better explore these trends, a large sample of carefully 
selected globules has to be studied. 

It is important to mention that we also used a slightly different fitting procedure, in 
which the theoretical profile was convolved with a Gaussian kernel of FWHM = 20", the 
resolution of the extinction maps. We obtained that, while the temperature remains nearly 
the same, the stability parameter ^^nax seems to be greater in all cases. Nevertheless, their 
associated uncertainty is also greater for all fittings carried out with kernel convolution, and 
in more than 60% of the cases the error bars overlap. Also, in all the fittings the reduced 
test is always bigger for convolved profiles, suggesting that there is no need to convolve 
a theoretical Bonnor-Ebert profile with a Gaussian kernel for fitting the dense cores of our 
sample using 2MASS data. 



5.4. Modelling Results and Discussion 



Figure [8] shows the stability parameter ^max versus the center-to-edge density contrast 
Pel Pr for the globules of our sample. Gray data correspond to starless globules and black 
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Table 6. Physical Parameters of Dense Cores in Bok Globules from Bonnor-Ebert Fitting 



Name 


^max 


T 


Tic 


M 


Pext/ks 










(K) 


(10" cm-3) 




(10* cm-3 K) 


(mag) 






Starless Globules 


016 A 


5.6 ± 0.1 


24.4 ± 1.4 


18.8 


3.1 


46.0 


18.2 


0.90 


016 B 


8.2 ± 0.2 


8.5 ± 1.2 


16.9 


1.1 


5.6 


10.8 


0.62 


044 A 


8.1 ± 0.5 


8.3 ± 1.0 


11.2 


1.2 


3.7 


8.7 


0.45 


044 B 


3.9 ± 0.7 


10.4 ± 1.7 


5.1 


1.0 


11.6 


5.6 


0.95 


053 


4.7 ± 0.5 


16.6 ± 1.5 


7.5 


2.1 


18.7 


9.1 


0.61 


059 


5.3 ± 0.4 


21.7 ± 1.2 


19.6 


2.4 


48.6 


17.3 


1.84 


074 


6.7 ± 0.5 


9.3 ± 1.2 


8.7 


1.3 


5.3 


7.9 


0.30 


075 


6.5 ± 1.6 


8.5 ± 1.3 


4.5 


1.6 


2.6 


5.4 


0.08 


111 


4.9 ± 0.2 


25.2 ± 1.8 


5.5 


4.9 


19.0 


9.7 


0.86 


113 


6.1 ± 2.0 


9.3 ± 2.0 


7.1 


1.3 


5.4 


7.0 


0.28 


133 


8.2 ± 0.8 


11.0 ± 3.0 


5.1 


2.8 


2.2 


6.8 


0.06 


144 A 


4.3 ± 0.1 


19.4 ± 1.6 


7.7 


2.3 


27.2 


9.7 


0.53 


144 B 


3.4 ± 0.3 


20.7 ± 1.7 


6.0 


1.9 


35.6 


8.1 


0.58 


145 


5.1 ± 0.3 


9.7 ± 1.1 


5.7 


1.2 


6.9 


6.2 


0.24 


IRAS Globules 


034 


4.5 ± 0.7 


7.9 ± 1.0 


5.1 


0.8 


6.6 


5.1 


0.72 


058 


3.8 ± 0.2 


22.7 ± 1.3 


12.5 


1.9 


65.8 


12.8 


3.02 


117 A 


3.9 ± 0.4 


8.9 ± 1.3 


4.8 


0.8 


9.4 


5.0 


0.27 


117 B 


3.8 ± 0.3 


7.1 ± 1.6 


2.7 


0.7 


4.5 


3.4 


0.21 


121 


9.4 ± 0.6 


5.1 ± 0.8 


13.3 


0.6 


1.9 


7.6 


0.46 


126 


6.8 d= 0.7 


10.1 ± 1.9 


4.2 


2.2 


2.7 


5.8 


0.08 


138 


6.4 ± 0.4 


10.7 ± 1.5 


4.9 


2.1 


3.8 


6.3 


0.33 


139 


6.4 ± 0.5 


17.6 ± 1.3 


15.6 


2.5 


20.0 


14.4 


1.95 


140 A 


4.8 ± 0.7 


12.5 ± 1.4 


7.3 


1.5 


13.0 


7.8 


0.49 


140 B 


7.0 ± 1.7 


7.2 ± 1.1 


10.7 


0.9 


4.5 


7.7 


0.44 


140 C 


6.1 ± 0.9 


10.7 ± 1.1 


12.0 


1.3 


10.5 


9.8 


0.97 


148 


4.0 ± 0.4 


10.6 ± 1.8 


3.5 


1.2 


7.8 


4.7 


0.16 


149 A 


5.1 ± 0.3 


12.7 d= 1.2 


10.2 


1.4 


16.1 


9.5 


0.75 


149 B 


6.3 ± 0.6 


6.7 ± 1.0 


8.2 


0.8 


4.1 


6.4 


0.24 



Table 7. Statistics for Dense Cores in Bok Globules 



All Globules Starless Globules IRAS Globules 



Physical Parameter 


Mean 


SD 


Median 


LQ 


UQ 


Mean 


SD 


Median 


LQ 


UQ 


Mean 


SD 


Median 


LQ 


UQ 
U 


R (pc) 


0.07 


0.02 


0.07 


0.06 


0.07 


0.08 


0.02 


0.07 


0.07 


0.07 


0.07 


0.02 


0.06 


0.06 


0.07^:5 


^max 


5.7 


1.6 


5.5 


4.5 


6.6 


5.8 


1.6 


5.5 


4.8 


6.7 


5.6 


1.6 


5.6 


4.1 


6.4 1 


r(K) 


12.6 


5.9 


10.5 


8.5 


16.9 


14.5 


6.5 


10.7 


9.3 


20.4 


10.8 


4.7 


10.4 


7.4 


12.1 


ric (10* cm-3) 


8.7 


4.7 


7.4 


5.1 


11.4 


9.2 


5.3 


7.3 


5.6 


10.6 


8.2 


4.2 


7.8 


4.8 


11.7 


M (Mq) 


1.7 


0.9 


1.4 


1.0 


2.1 


2.0 


1.1 


1.7 


1.3 


2.3 


1.3 


0.6 


1.2 


0.8 


1.8 


Pe:ct/kB (lO* Cm-3 K) 


14.6 


16.2 


7.4 


4.4 


18.8 


17.0 


16.2 


9.3 


5.3 


25.2 


12.2 


16.3 


7.2 


4.2 


12.4 


(mag) 


8.5 


3.6 


7.8 


6.1 


9.7 


9.3 


3.9 


8.4 


6.8 


9.7 


7.6 


3.1 


7.0 


5.3 


9.0 
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data to IRAS globules. This figure shows that there is no difference between starless and 
IRAS globules with respect to the stability parameter. Taking into account the uncertainties 
in imax^ we can see three groups in this plot: a first group showing clearly stable states, with 
imax = 4.5 ± 0.7 (8 starless + 7 IRAS), a second group distributed around the critical state 
with ^rnax = 6.5±0.2 (3 starless + 6 IRAS), and a third group showing clearly unstable states, 
with ^rnax = 8.5 ± 0.6 (3 starless + 1 IRAS). Remarkably, the first group of stable globules 
has a temperature of T = 15 ±6 K, and the critical plus unstable groups have T = 10 ±3 K. 
This tendency of decreasing temperature with increasing stability parameter can also be 
seen in Figure [9|, where both parameters are plotted with their respective uncertainties. The 
same behavior is shown by the external pressure. 



These results show a difference with respect to the sample of lKandori et al.l (120051 ): they 
found that more than half of their starless globules are located near the critical state, which 
is not our case. Secondly, their Bonnor-Ebert temperature distribution (not shown by these 
authors) against stability parameter did not show any tendency of decreasing temperature 
with increasing ^max, as in our sample. 

Figure [TD] shows histograms of the logarithmic density contrast, with the vertical dashed 
line corresponding to the critical Bonnor-Ebert sphere. The left histogram corresponds to 
the whole sample, showing that there is a large number of cores in the stable regime, followed 
by an abrupt decrease of cores in the unstable regime. The same behavior is seen for starless 
globules (upper-right) and for IRAS globules (lower-right), but the decrease in the number of 
cores for IRAS globules is steeper than for starless globules. This seems to be in agreement 
with the fact that globules with density contrasts larger than 14 must have shorter lifetimes. 
The opposite happens when the density contrast is smaller than the critical value. In order 
to state whether this is a real trend or not a large sample of globules has to be considered. 
The peak of the starless globules occurs just below the critical value, while for IRAS globules 
th e peak is aroun d the c ritical value. This is in contrast to the globule distribution showed 



by iKandori et al.l (l2005l . see their Figure 6), where the majority of starless globules are in 



the unstable regime and the star-forming globules are well in the unstable regime. 



Similarly to lKandori et al.l (120051 ). we plot in Figure [TT] the correlations between density 
contrast and the different parameters of the globules derived from the Bonnor-Ebert fitting 
in a logarithmic scale. The dotted lines correspond to the variation of each parameter for 
Bonnor-Ebert spheres with minimum and maximum values of T and Pext {Tmin = 5.1 K and 
Pext,min = 1-9 X 10"^ cm~^ K; Tmax = 25.2 K and Pext,max = 65.8 X 10'^ cm~^ K), as indicated 
in Figures [TT] (a) and (b) for the temperature and external pressure, respectively. With these 
two extreme Bonnor-Ebert spheres, we see in Figures [TT] (c)-(f) that we can cover nearly all 
the points in the plots, which demonstrate that the dispersions in the distributions of the 



- 23 - 



parameters can be accounted by two Bonnor-Ebert spheres, one with fixed minimum values 
of T and Pexti and other with fixed maximum values. 

Finally, if the IRAS globules are star-forming globul es, then our sample of Bok globules 
do not show an evolutionary sequence, as is the case of iKandori et al.l (120051 ). where their 
sample of globules suggest an evolutionary sequence in the sense that the starless globules are 
distributed around a critical Bonnor-Ebert sphere and the star-forming globules are highly 
unstable. 



6. The Infrared Sources 

To study the properties the of IRAS sources associated with the globules of our sample, 
we looked for fiux measurements in the following catalogs: 

• DENI^Ifor band I (0.79 /xm) 

• 2MASS for bands J (1.24 /xm), H (1.66 /im) and K (2.16 /im) 

• MSxIfor bands A (8.28 /im), C (12.13 /im), D (14.65 /im) and E (21.34 /im) 

• IRAS for bands 12 /im, 25 /im, 60 /im and 100 /im 

To determine the nature of the IRAS sources embedded in the globules, we used two 
criteria to segregate the spectral classes: in §6.11 we build the spectral energy distribution 
(SED) for the IRAS sources and calculated the spectral index, and in §6.21 we calculated the 
bolometric luminosity and temperature for the sources. 



6.1. Spectral Energy Distribution 

The shape of the broadband spectrum of a YSO depends on the nature and distribution 
of the circumstellar material. This implies a correlation between the SED and the evolution- 
ary state of the YSO, where the initial protostellar stages, during which a stellar embryo is 
surrounded by large amounts of collapsing circumstellar material, have very different infrared 



''Deep Near- Infrared Survey ( Epchtein et al.lll99!^ 
^Midcourse Space Experiment ( Price et al.ll2001 ) 
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signatures than the more advances stages of pre-main sequence, where most of the original 
material has been incorporated to the young star (see lLadalll999l ). 



Class (lAndre et al.lll993l ) and I fj Adams et al.lll987l : lLadalll987l ) sources are character- 
ized by SEDs with peaks in the far-infrared and sub-milimeter, as a result of the emission 
of a cold dust collapsing envelope. These objects are the most deeply embedded and ob- 
scured. Once the collapse of the envelope ceases, the dusty disk can still produce substantial 
emission in the infrared. This part of the evolution gives raise to Class II and III sources 
( lAdams et al.l 119871 : iLadal 119871 ). whose SEDs are characterized by peaks in the opitcal and 
near-infrared, where stellar fotospheric emission dominates. These objects correspond to T 
Tauri stars (CTTS and WTTS). 



Based on the shapes of SEDs of observed infrared sources, I Adams et al.l (119871 ) defined 
the spectral index a between 2-25 /xm as: 



a 



dlogjXF,) 
rflog(A) 



(27) 



where F\ is the flux corresponding to wavelength A. The different spectral indices computed 
in regions of star formation are interpreted as an evolutionary sequence, obtaining the fol- 
lowing classification: a > for Class I sources, > a > —2 for Class II sources, and a < 
—2 for Class III sources. 



The SEDs of the IRAS globules are shown in Figure [121 To calculate the spectral index 
we fitted a straight line using the last squares method to the points between 2.2 and 25 fim, 
with the requirement that these points have a determined uncertainty. Fluxes that are upper 
limits (shown as arrows in Figure [T^ where not taken into account. Table [H] shows the name 
of the globule in column 1, the associated IRAS sources in column 2, the index a in column 
3 and the corresponding spectral classification according to that index in column 4. It was 
possible to determine the spectral index for 8 of 13 IRAS sources, of which 6 correspond to 
Class I, one to Class II and one to Class III. 



6.2. The BLT Diagram 

A useful tool for studying the evolutionary state of a young object is th e graph of the 
bolometric luminosity (Lboi) vs. bolometric temperature (Tboi) (BLT diagram). iMyers fc Ladd 
(119931 ) defined and used the BLT diagram to study the evolutionary state of YSOs in the 
Taurus-Auriga complex. T^oi is defined as the temperature of a blackbody which has the 
same mean frequency of the observed continuum spectrum. A main-sequence star whose 
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Table 8. Spectral Classification of YSOs in Bok Globules 



Name 


IRAS 


a 


Class^ 


Lbol 


Tbol 












(Le) 


(K) 




034 


08250-5030 






1.0 


29 





058 


10471-6206 


1.3 


I 


2.2 


32 





117 


16029-4548 


-1.0 


II 


4.2 


742 


II 


121 


16549-5030 






0.6 


227 


I 




16554-5031 


1.7 


I 


1.3 


98 


I 


126 


16009-3927 


-2.5 


III 


37.6 


2478 


II 


138 


17159-4324 






3.3 


28 





139 


17172-4316 


0.7 


I 


2.8 


123 


I 




17169-4314 






1.9 


90 


I 


140 


17193-4319 


1.5 


I 


1.9 


35 







17195-4320 






15.4 


1676 


II 


148 


17011-3613 


1.6 


I 


1.6 


31 





149 


17012-3603 


1.1 


I 


1.5 


39 






'^Spectral classification according to index a 
'^Spectral classification according to Tboi 
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spectrum is a blackbody has Tboi = Teff. A pre-main sequence star with infrared excess has 
a wider and redder spectrum than a blackbody, and therefore T^oi < T^ff. An embedded 
infrared source is very reddened to have an optical measur e of T^/ /, its spec trum is much 
wider than a blackbody, and its T^oi is very small (< 500 K, IChen et al.lll995l ). 



To derive the bolometric luminosity for the IRAS sources, we integrated their SEDs 
across the whole range in wavelengths: 



^bol 



XFx dlogX. 



(28) 



In most cases, the lowest available wavelength to make the integration is the J band. As the 
SEDs for all the sources grow up as the wavelength increases (see Figure [T2l) . we ignored the 
terr a between A = and the J band, which is very small compared to the total luminosity 
(see 



Baba et al.ll2006l ). Therefore, the integration was divided into two parts: 



^bol 



9.27TD' 



A2 poo 

XFxdlogX+ / XFxdlogX 

Ai J X2 



(29) 



where Ai = J and A2 = 100 /im, the longest available wavelength. To compute the second 
part of the integral, as we do not know the shape of the spectrum be tween A2 and infinity, 
we assumed that the SED is proportional to A~^ ( IWilking et al.lll989l ). 



The bolometric temperature is defined by iMyers fc Laddl (119931 ) as: 



hoi 



C(4) hV 
4C(5) k 



B 



1.25 X 10"^^I7(KHz"^), 



(30) 



where Ci"^) is the Riemann zeta function of argument n and the mean frequency v is the 
ratio of the first and zeroth moments of the source spectrum: 



V 



uFiy dv 



(31) 



Once we calculated these parameters, we constructed the BLT diagram. IChen et al. 



( I1995I ) established a correlation between the temperature ranges in the BLT diagram and 
the evolutionary class of a YSO, getting as a result the following classification: T^oi < 70 K 
for Class 0, 70 K < T^oi < 650 K for Class I, 650 K < < 2800 K for Class II, and T^oi > 
2800 K for Class III sources. The values of Tboi, L^oi and spectral classification according 
to this criterion are shown in columns 5, 6 and 7 of Table [HI respectively. Out of 13 IRAS 
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sources, 6 were classified as Class 0, 4 as Class I, and 3 as Class II . The BLT diagram is 



shown in Figure [T3l where the solid line is the ZAMS from lSiess et al.l (120001 ) and the vertical 



dashed lines represent the intervals corresponding to different spectral classes. 

The mean value of the Bonnor-Ebert stability parameter is ~ 5 for each spectral class, 
suggesting that there is no difference between the average value of the stability parameters 
of the cores and the spectral class of their associated IRAS sources. 



7. BHR 138: A Collapsing Globule? 

In this section, we will use the millimetric observations data of the ^^CO and C^^O 
molecules, in the rotational transition J = 1 ^ 0, to study BHR 138 and BHR 149 (core 
A). B oth globules are part of a survey of more than 30 condensations identified in Scor- 



pius ( IVilas-Boas et al.ll2000l ). Among all these condensations, only BHR 138 presents a 
blue- asymmetric profile, while BHR 149 shows a typical profile observed in the others con- 
densations. The presence of a blue-asymmetric profile in the spectrum of BHR 138 provides 
evidence of an eventual gravitational collapse. No spectral line data were available to study 
the other globules. 



7.1. Reduction of the Spectra 

The reduction of ^^CO and C^^O spectra was made with the spectral lines reduction 
program Drawspec§|. The lines were adjusted with gaussians. In the case of BHR 138, the 
fitting was made with two independent gaussians. For m ore details about the re duction of 
the spectra and the derivation of the line parameters, see IVilas-Boas et al. J1994I ). 



Figure [H shows the ^^CO and C^^O spectra for BHR 138 and BHR 149. The histogram 
corresponds to the observed spectrum and the grey gaussians are the fitting made with 
Drawspec. For the ^'^CO line in BHR 138, the resulting gaussian is the sum of the two thin 
gaussians. Table [9] shows the parameters extracted from the line fitting. Column 1 is the 
antenna temperature corrected from atmospheric attenuation, column 2 is the velocity of 
the line center, and column 3 is the FWHM of the line. 



http: //www. cv.nrao . edu/'^{}hliszt 
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7.2. Modelling of the ^^CO Line 



Comparing the spectra of ^'^CO for BHR 138 and BHR 149 in Figures [T4l (a) and (b), 
respectively, we can see clearly the asymmetry in the line of BHR 138, while BHR 149 shows 
a perfectl y symrnetric line. The spectrum of BHR 138 does not presents two distinct peaks 



[see, e.g.. 



Evans 1999[). but do presents a typical profile of collapsing clouds denominated 



red-shoulders" by lMyers et al.l (Il996l ). 



Given the blue-asymmetric profile shown by the ^'^CO line in BHR 138, one of the 
fundamental questions is to know if this is really a profile indicative of collapse. One of 
the possibilities to explain that profile will be the superposition of distinct clouds in the 
line-of-sight or rotational effects of the clouds. Analysis of the field of the globule (Figure [2]) 
shows a completely isolated object and the antenna beam contains only the globule. Similar 
situation occurs with BHR 149, which is also isolated and the profile does not suggests any 
asymmetry in the lines. Analysis of 30 globules in Scorpius also does not show profi les with 
asymmetries of the kind identified in BHR 138. According to lVilas-Boas et al.l (120001 ). only 4 
condensations of their sample observed in that region present contamination for superposition 
of clouds in the antenna beam. Generally, the components that contaminate the CO emission 
in that region have very wide line widths (> 4 km s~^) at different radial velocities. 

Because the ^^CO and C^^O emissions are orig inated in the more external regions of the 
globules, with densities lesser than 10^ cm~^ (e.g.. lAlves et al.lll999l ). we would also expect 



to identify that asymmetry in the C^^O line. Nonetheless, the signal-to-noise ratio for the 
spectrum of that molecule does not allow to identify a red-shifted component with sufficient 
accuracy. More greater integrations will be necessaries to improve the signal-to-noise ratio. 



The first model, proposed by lMyers et al.l (Il996l ). is called two-layer and consists of two 
collapsing gas layers, one approaching the other with a velocity of twice the infall velocity 
Vin- Each layer has a constant excitation temperature: the front layer, moving away from 
the observer, has a temperature Tf equal to the background temperature (Tj = 2.7 K), while 
the rear layer, which is hotter, has a temperature T^. Assuming that both layers have the 
same velocity dispersion a and optical depth tq at the line center, iDe Vries fc MyersI (120051 ) 
obtained the brightness temperature emerging from the spectral line: 



ATb = J(T,)(1 - e-^^)e-^f + J{Tf){e-^r 



(32) 



where J{T) = To/[exp(To/T) — 1], Tq = hu/ks, is the transition frequency , and Tf and 
are the optical depth of the radiation coming from the front and rear layers, respectively: 
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Tf{v) = tq exp [-{v - vlsr - Vinf/2a'^] 
Tr{v) = To exp [~{v - Vlsr + vinf 120"^] 



(33) 
(34) 



where vlsr is the hne-of-sight velocity, assumed to be the same for both regions. Therefore, 
this model has five free parameters: 

TO,Vin,Tr.,a, VlSR- 



The second model, introduced in iDe Vries fc Myerd (l2005l ). called hill model, consists of a 
dense core with a peak excitation temperature Tp at the center, and a temperature Tq at 
the core edges. The J(T) profile decreases linearly from J{Tp) at the center to J(To) at the 
edges. The optical depth is tc and the infall velocity is vc- Assuming To = 2.7 K, this model 
has five free parameters: 



'rc,vc,Tp,a, vlsr- 



For a graphical representation and different variants of the models of iDe Vries &: Myers 
(120051 ). see their Figure 1 and their Tables 1 and 2. 



7.3. Modelling Results 

To perform the fitting of the parameters for each model, we used a routine developed 
by De Vries fc Myers ( 2005 ). obtained from the site of C. H. De Vrie^. Starting values must 
be supplied to the program to find the best fitting parameters. Table [10] shows the results 
obtained with the two-layer model, and Table [TT] shows the results corresponding to the hill 
model. The model that best fit the observational data for BHR 138, according to the 
value, is the two-layer model. Figure [15] shows the fitting obtained with this model, where 
the grey curve represents the theoretical fit, which suggests a velocity collapse of 0.25 km s~^. 
Considering that rotation was not taken into account, the derived collapse velocity should 
be seen as upper limit. 



Possible rotational effects of the globules (iBarranco fc Goodman! 119981 ) would be eval- 
uated from maps of optically thin molecular transitions (e.g., N2H"'"). Studies of rotation 



" jhttp: //cf a- www. harvard. edu/'^{}cdevries/analyticJ.nf all .html 
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Table 9. Observed Parameters of Spectral Lines 



n Vlsr Ay 

(K) (km s"^) (km s"^) 



BHR 138 



CO 2.57 ± 0.11 -7.62 ± 0.01 0.78 ± 0.03 
1.32 ± 0.10 -7.27 ± 0.04 2.18 ± 0.09 



Ci^O 


0.46 ± 0.05 


-7.55 ± 0.05 


0.82 ± 0.11 


BHR 149 


"CO 


3.12 ± 0.03 


1.11 ± 0.01 


1.02 ± 0.01 


C^^O 


0.31 ± 0.03 


1.08 ± 0.05 


1.01 ± 0.12 



Table 10. Two-Layer Model Fitting for BHR 138 



To Vin 


Tr CF VLSR 




(km s~^) 


(K) (km s-i) (km s'^) 




1.6 0.25 


10.9 0.6 -6.9 


1.2 



Table 11. Hill Model Fitting for BHR 138 



Tc vc Tp a VLSR 

(km s-i) (K) (km s-^) (km s-^) 



4.2 0.49 7.7 0.3 -7.1 3.8 
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in starless globules through the transition N H.s (1,1) show maximum velocity gradients, du e 



to rotation, of the order of 1.4 km s ^ pc ^ (IBarranco &: Goodmarull998l: ISwift et al 



2005), 



2002) 



and a survey of 12 starless globules and 14 protostellar globules in N2H"'" (jCaselli et al. 
suggest typical velocity gradients of 2 km s~^ pc~^ in both kind of globules. 

Adopting a mean velocity gradient due to rotation of 1.4 km s~^ pc~^ as being typical 
for Bok globules, it is possible that the collapse velocity derived for BHR 138 could be 
overestimated by 40% due to rotational effects that were not corrected in this work. Of the 
factors that would affect this conclusion, rotation seems to be the most relevant component. 
Taking that gradient and the size of BHR 138, rotation could be contaminating the spectrum 
with 0.14 km s~^, assuming the angular momentum perpendicular to the line-of-sight. In 
that case, the collapse velocity would be ~ 0.1 km s~^. 



In a sample of 53 starless cores observed in N2H"'", iLee et al.l (1200 ll ) fou nd collapse 
velocities of ~ 0.1 km s~^ using the two-layer model. Using observations of CS. ISwift et al. 



(120061 ) obta ined a collapse velo city of 0.15 km s~^ for the starless core L1551-MC with the 
hill model. iMyers et al.l (119961 ) . utilizing the two-layer model, derived a collapse velocity of 
0.35 km s^^ for the dense core L1251B, which harbors a protostar. 



8. Summary 

In this work, we studied a sample of 21 isolated Southern Bok globules selected among 
169 globules catalogued by BHR. In this sampl e, 11 are starless an d 10 have associated IRAS 



point sources. Employing the NICE method of iLada et al.l (119941 ). we used the 2MASS PSC 
to construct visual extinction maps for the Bok globules with 20" resolution. This enabled 
us to detect dense cores and substructures embedded in the globules, which were modelled 
with Bonnor-Ebert spheres to derive their physical parameters and investigate their stability. 
Two methods were used to determine the spectral class of the IRAS sources associated to the 
globules. In addition, we searched for spectral lines observations data obtained toward these 
globules to explore whether they have or not evidence of collapse. We found observations 
of ^^CO and C^^O J = 1-0 lines toward two globules of the sample showing evidences of 
collapse in one of them. Distances were estimated to all these globules from plots of color 
excess versus distance. The main results of this paper are summarized as follows: 

1. Dense cores were identified embedded in the globules, with mean volumetric densities 
greater than 10^ cm~^ and masses between ~ 1-4 M©. 



2. We inferred distances to several globules different from that ones obtained by BHR in 
spite of using the same technique. Due to the uncertain nature of this method, in order 
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to improve the distance determination, we used smaller areas and a catalog containing 
information for approximately 10 times more stars than the catalog utilized by BHR. 

3. We found that there is a large number of cores in the stable regime, followed by an 
abrupt decrease of cores in the unstable regime, but this decrease is steeper for IRAS 
globules than for starless. The globules that are in the stable state have a temperature 
oi T — 15 ± 6 K, while the globules that are in the critical and unstable states have a 
remarkably different temperature ofT'=10±3K. 

4. We found that, according to the bolometric temperature, 6 IRAS sources are classified 
as Class 0, 4 as Class I, and 3 as Class II. The spectral index also produced similar 
results. No correlation was identified between the Bonnor-Ebert stability parameter of 
the condensations and the spectral class of the associated IRAS sources. 

5. By means of the ^'^CO J = 1-0 molecular line, we identified a blue-assymetric line profile 
toward BHR 138. Using the two-layer collapse model, it was possible to estimate an 
infall speed of 0.25 km s~^. However, taking into account possible rotational effects, 
this value could be as small as ~ 0.1 km s~^. 



We are grateful to an anonymous referee for useful comments and suggestions that 
improved the manuscript. G. A. R. is grateful to a PCI-INPE/CNPq fellowship and a 
CAPES grant. This publication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts and the Infrared Processing 
and Analysis Center, funded by the National Aeronautics and Space Administration and the 
National Science Foundation. 
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Fig. 1. — Extinction map for B68 constructed using 2MASS data with resolution 20". Con- 
tours start a t Ay = 4 mag and increases in steps of 2 to 26 mag. Compare with the map of 



Alves et al.l (120011 ). 
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Fig. 2. — Visual extinction maps constructed with 2MASS data using the NICE method. 
Each map has ~ 15' x 15' and resolution h = 20". 
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Fig. 2. — continued. 
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Fig. 3. — Visual extinction maps constructed whith 2MASS data and the star counts method. 
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Fig. 4. — Galactic distribution of the Bok globules of our sample. The numbres are the 
names of the globules given in Table [Tj and the circles represent the 3° regions used to 
search for stars in the SKY2000 catalog. In some cases (58-59, 74-75, 138-140, and 148-149), 
the same region was used. 




Fig. 5. — Color excess versus distance for stars in the neighborhood of each globule. The 
dashed vertical lines indicate the distances adopted in this work. In the upper left corner 
are given the distances estimated to the globules. 
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Barnard 68 




^mnx = 6.9 ± 0.3 

T = 16.8 ± 0.9 K 







20 40 60 80 100 

Radius (arcsec) 



120 140 



Fig. 6. — Azimuthally- averaged radial extinction profile for Barnard 68. The solid curve 
represents the Bonnor-Ebert theoretical profile fitting. 




Fig. 7. — Azimuthally-averaged radial extinction profiles for dense cores in Bok globules. The 
solid curve superimposed on each observed profile represents the theoretical Bonnor-Ebert 
profile fitting. 




Fig. 7. — continued. 
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Fig. 8. — Relation between ^max and the center-to-edge density contrast for dense cores in 
Bok globules. Grey dots correspond to starless globules and black dots to IRAS globules. 
The vertical dashed line denotes the critical value ^^ax = 6.5. 
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Fig. 9. — Relation between stability parameter and temperature for dense cores in Bok 
globules. The stable globules have T — 15 ± 6 K and the critical plus unstable ones have 
T = 10 ± 3 K. Grey dots correspond to starless globules and black dots to IRAS globules. 
The vertical dashed line denotes the critical value ^max — 6.5. 
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Fig. 10. — Histograms of the logarithmic density contrast for the whole sample of Bok 
globules (a), for starless globules (b), and for IRAS globules (c). The vertical dashed line 
denotes the critical value of the density contrast, pd Pr =14. 
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Fig. 11. — Correlations between density contrast and tlie parameters determined from tlie 
Bonnor-Ebert fitting of dense cores in Bok globules, (a)-(f) Plots for temperature, external 
pressure, central extinction, central density, radius and mass, respectively. 




Fig. 12. — Spectral energy distributions for infrared IRAS sources associated to Bok globules. 



-47- 




Fig. 13. — BLT diagr am of infrared IR AS sources associated to Bok globules. The solid line 
indicates the ZAMS (ISiess et al.ll2000[). Th e vertical dashed lines determine the limits of the 
different spectral classes ( Chen et al.l 1995 ). 
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Fig. 14.— i^CO and C^^O line profiles for (a) BHR 138 and (b) BHR 149. 



-49- 



5 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




' ' 

-11 -10 -9 -8 -7 -6 -5 -4 

Vlsr (km/s) 



Fig. 15.— Modelling of the ^^CO line profile for BHR 138. The histo gram represents the 
observational data and the grey curve corresponds to the two-layer fit model, which provides 
an infall velocity of 0.25 km s~^. 
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